Optical activity of noncentrosymmetric metals 
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We describe the phenomenon of optical activity of noncentrosymmetric metals in their normal 
and superconducting states. The found conductivity tensor contains the linear in wave vector off 
diagonal part responsible for the natural optical activity. Its value is expressed through the ratio 
of light frequency to the band splitting due to the spin-orbit interaction. The Kerr rotation of 
polarization of light reflected from the metal surface is calculated. 



PACS numbers: 78.20.Ek, 74.25. Nf, 74.20.Fg 

I. INTRODUCTION 

The metals without inversion symmetry have re- 
cently become a subject of considerable interest mostly 
due to the discovery of superconductivity in CePtsSi. 1 
Now the list of noncentrosymmetric superconductors 
has grown to include Ulrl, CeRhSi-A CelrSi^, Y 2 C 3 5 , 
Li2(Pdi_ 2; ,Pt a ;)3B 6 , KC^Oe 7 -, and other compounds. 
The spin-orbit coupling of electrons in noncentrosymmet- 
ric crystal lifts the spin degeneracy of the electron energy 
band causing a noticeable band splitting. The Fermi sur- 
face splitting can be observed by the de Haas- van Alphen 
effect discussed theoretically in the paper Ref. 0. The 
band splitting reveals itself in the large residual value of 
the spin susceptibility of noncentrosymmetric supercon- 
ductors at zero temperature^ It also makes possible the 
existence of nonuniform superconducting states that can 
be traced to the Lifshitz invariants in the free energy*^ 
Another significant manifestation of the band splitting is 
the natural optical activity. 

The natural optical activity or natural gyrotropy is 
well known phenomenon typical for the bodies having no 
centre of symmetry-ii. The optical properties of a nat- 
urally active body resemble those of the magneto-active 
media having no time reversal symmetry. It exhibits dou- 
ble circular refraction, the Faraday and the Kerr effects. 
In this case, the tensor of dielectric permeability has lin- 
ear terms in the expansion in powers of wave vector 



/(<*>, q) 



(1) 



where 7^ (w) is an antisymmetric third rank tensor called 
the tensor of gyrotropy. 

The spacial dispersion term in permeability has been 
derived by Arfi and Gor'kovii 2 - in frame of model where " a 
conductor lacking a center of inversion is simulated by an 
ordered arrangement of impurities whose scattering po- 
tential is asymmetric" (seealsoi^). We shall be interested 
in gyrotropy properties of a clean noncentrosymmetric 
metal. In metals, it is more natural to describe optical 
properties in terms of spacial dispersion of conductivity 
tensor having the following form: 



Here, the gyrotropy tensor Ayi(tj) is an odd function of 
frequency. 

The gyrotropic tensor has the most simple form in the 
metals with cubic symmetry. In this case, the Drude 
part of the conductivity tensor is isotropic tjij(ui,0) — 
a(uj)Sij and the gyrotropic conductivity tensor A^fcz (<^) = 
\{w)eiki is determined by the single complex function 
\(ui) = \'(uj) + such that a normal state density 

of current is 



e 9E „ , „ 
J = — Trr + ctE + ArotE. 

in at 



(3) 



In time representation A is an operator being an odd 
function of operation of time derivative d/dt. 

In the superconducting state besides ArotE the gy- 
rotropic part of current density contains also an addi- 
tional term proportional to the magnetic induction B 



j„ = ArotE + vB. 



(4) 



;(w,q) 



/(w,0) 



i\iji(uj)qi. 



(2) 



Here v = u(T) is a constant coefficient being equal to 
zero in the normal statei 14 ' 15 

In this paper we present the derivation of current re- 
sponse to the electro-magnetic field with finite frequency 
and wave vector valid for the normal and the supercon- 
ducting state of noncentrosymmetric metals with cubic 
symmetry. We find that the gyrotropy conductivity A is 
directly proportional to the ratio of the light frequency 
to the energy of the band splitting. Then making use 
the Maxwell equations and corresponding boundary con- 
ditions at the surface of noncentrosymmetric metal de- 
rived in the Section V we calculate the Kerr rotation of 
polarization of light reflected from the surface of metal. 



II. CURRENT RESPONSE TO 
ELECTRO-MAGNETIC FIELD 

The current response of a clean metal to the electro- 
magnetic field at finite q and ui can be written following 
the textbook procedure^. In application to our situation 
one has to remember that due to spin-orbital coupling 
determined by the dot product of the Pauli matrix vec- 
tor (7 and pseudovector 'y(k), which is odd in respect to 



2 



momentum 7(— k) = — 7(k) and specific for each non- 
centrosymmetric crystal structur e) 9 ' 17 , all the values such 
as single electron energy 



where 



&*/3(k) = (e(k) - n)5 a p + 7(k)<T a/3 , 



(5) 



velocity v Q/ 3(k) = 9£ Q/ 3 (k) /<9k, the inverse effective 
mass (m~ 1 ) Q ^ = d 2 t; a fi(k) / 'dkidkj, the Green func- 
tions G , ct ^(Ti,k;r 2 ,k / ) = -(T T a ka (Ti)aj c , /3 (r 2 )) and 
F Q( g(ri,k;T 2 ,k') = (T T a kQ (Ti)a_ k > ^(r 2 )) are matrices in 
the spin space. Taking this in mind, we obtain 



ji(u n ,q) = -—Tr [rh^fie + 



/ 7^4 T E {^(k)G (0) (^ + )%(k)G(°)(^ 



m— — 00 



(27ra) 3 

+^(k)F( )(^ + )^(-k)^+(°)(^_)}j A>„,q). (6) 

The transposed matrix of velocity is determined as 
v'(— k) = <?£*(— k)/9k. The arguments of the 
zero field Green functions are denoted as K± = 
(£} m ± oj„/2, k ± q/2). The Matsubara frequencies take 
the values VL m = 7r(2m + 1 — n)T and u n = 2imT '. 

One can pass from the spin to the band representation, 
where the one-particle Hamiltonian 

ff o = ^W( k ) a L fl k/3= E ^(k)4 A c kA (7) 

k k,A=± 

has diagonal form. Here, the band energies are 

e A (k)=e(k)- A1 + A|7(k)|, (8) 

such that two Fermi surfaces are determined by equations 
£ A (k) — 0. The difference of the band energies 2|7(k^)| 
characterizes the intensity of the spin-orbital coupling. 
The Fermi momentum with 7 = is determined by the 
equation e(kp) = ef- 

The diagonalization is made by the following transfor- 
mation 



with the coefficients 



(9) 



A=± 




V2|7l(l7l+A 7z )' 



forming a unitary matrix ii(k). 

The zero field Green functions in the band representa- 
tion are diagonal and have the following form: 9 



G AA ' ( w « , k) = S\\'G\ (uj n , k) , 

AA' 



F^,(uj n ,k)=5 X yF x (uj n ,k), 



(10) 



G A (w„,k) = - 
F x (u n ,k) = 



ihlUn + £ A 



and 



*A(k) = -A 



^ + ei + |A A (k)|2' 

t A (k)A A (k) 
7a; (k) -ily(k) 



(11) 



V / 7 2(k)+ 7 2(k)' 



The functions A A (k) are the gaps in the A-band quasipar- 
ticle spectrum in superconducting state. In the simplest 
model with BCS pairing interaction v g (k, k') = —V g , the 
gap functions are the same in both bands: A+(k) = 
A_(k) = A and we deal with pure singlet pairing 18 . 

Transforming the Green functions using eqn. ([§]) into 
the band representation^, we obtain 

Tr{v t (k)G (0 \K + )v 3 (k)G^(K_) 
+v l (k)F^ (K + )v t 3 (-k)F+^ (K_)} = 
W++ , i (k)G +U++J (k)G+ + v ++!i (k)F + v ++d (-k)Fl + 
v--,iOL)G-v--j(\L)G-+v--,i{k)F-v- + j(-k)Fl + 
w+_ )i (k)G_t;_ + ^(k)G + + u+_, i (k)F_v_ +ii (-k)4 + 
v_+, i (k)G + « + _ ii (k)G_ + w_ +li (k)F + « + _ iJ -(-k)fl (12) 

as the trace of the matrices in eqn. ©. For the brevity, 
we omit here the arguments of the Green functions. They 
are the same as in the upper two lines. The matrix ve- 
locity in the band representation is 



v AV (±k) = u[ a (±k)v al3 (±k)u fjx ,(±k) = 

0&(±k) d 7l (±k) 

-o\\> H 7^- — Ti,\x'{±k), (13) 



9k 



9k 



where "f(k) = iv (k)o"fi(k) are hermitian matrices. We 
have neglected^ the difference between u(k) and u(k ± 

q/2). 

The explicit expressions for the r(k) matrices are 



% 


7x72+''7m 

7-L 


7X 


7v7*+*7s 
7-L 


7m 7 2 -*7x 

7-L 

-ly 


7x \ 




7 I 





(14) 



^ -7, 



where 7 = 7/I7I, 7 _l = y 7 2 +7y- All the diagonal ele- 
ments of these matrices are odd functions of the momen- 
tum components. Correspondigly, the products of the 
diagonal components of velocity matrices (|13[) are even 
functions of momentum. Hence, these terms in eqn. (|12[) 
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can produce only the terms proportional to even powers 
of the product kq being not responsible for the gyrotropy 
properties. 

The off diagonal elements of r(k) have the mixed par- 
ity. So, the dispersive terms proportional to the odd 
powers of the product kq can arise only from the part 
of eqn. (fT2|) containing the off diagonal elements of r(k) 
matrices. Hence, for the calculation of gyrotropy of con- 
ductivity, only the last two lines in eqn. (TT2|) consisting 
of interband terms are important. They are equal to 



dkj dk. 



7i,_ + r* ) _ + [G + G_-F + f!]}. 



(15) 



Using the identity 



where 7; = one can rewrite the above expression 

as 



dh dk. 



-{(Sim - ia m )[G-G+ + G+G_ - F_M (16) 



-F+Fl] + ie lmn %[G-G+ - G+G_ - F_F| + F+f!_}}. 

Starting this point we need the explicit form of spin- 
orbit coupling vector 7(k). Its momentum dependence 
is determined by the crystal symmetryj&il For the cu- 
bic group G = O, which describes the point symmetry 
of Li2(Pdi_ 2; ,Pt 2 ;)3B, the simplest form compatible with 
the symmetry requirements is 

7 (k) - 70 k, (17) 

where 70 is a constant. For the tetragonal group G = 
C4 V , which is relevant for CePtaSi, CeRhSi3 and CeIrSi3, 
the spin-orbit coupling is given by 

7(k) = 71 (kyx - k x y) + 7|| k x k v k z (k 2 x - k 2 )z. (18) 

The gyrotropy current j g , which is linear with respect 
to the wave vector q, originates from the last term in the 
eqn. (fl6|) . One can show that for the tetragonal crystal 
with the symmetry group G = C^, for the electric field 
lying in the basal plane the linear in the component of 
wave vector q part of conductivity is absent. In that 
follows we continue calculation for the metal with cubic 
symmetry where 7 = k sigwyo- We put 7 = k taking 
70 as a positive constant. Thus, we obtain for gyrotropy 



current* 2 ^ 



e 2 7 2 



(19) 



d 3 k 



■h 



(2tt?i) 3 

OO 

xT [G+(K + )G-(K-)-F + {K + )Fl{K_)- 

m— — 00 

-G_ (K + )G + (K_ ) + F_ (K + )Fl (K_ )] . (20) 

Let us find first the gyrotropy conductivity in the nor- 
mal state. 



III. GYROTROPY CONDUCTIVITY IN THE 
NORMAL STATE 

Substituting the Green function in the eqn. (|2"tJ| and 
performing summation over the Matsubara frequencies, 
we obtain 



d 3 k 
{2ttTi) : 



-h 



/(g_(k_))-/(g + (k + )) 
i&w„ + £-(k_)-£+(k + ) 

/(£ + (k_))-/(£_(k + )) 
iRw n + e+(k_)-f_(k+) 



(21) 



Here /(£±(k±)) is the Fermi distribution function and 
k± = k ± q/2. By changing the sign of momentum k — > 
— k in the first term under integral and making use that 
£a (k) is even function of k, we come to 



h 



d 3 k ^K+(k+) 



£_(k_)][/(£ + (k + ))-/(Mk_))] 



{2irh) 



(f+(k + )-e_(k_))a-(in Wn )2 

Analytical continuation of this expression from the dis- 
crete set of Matsubara frequencies into entire half-plane 
u) > is performed by the usual substitution iuj n — > 
uj + i/r. 

We shall work at frequencies smaller when the band 
splitting Tluj < "fo^F far from the resonance region Tiui ~ 
7o&f but still in the collisionless limit lot > 1 where one 
can decompose the integrand in powers of w 2 : 



I, =2 



d 3 k 



■fc*[/(£ + (k + ))-/(£_(k_))] 



+ 



(Hlo) 



€ + (k+)-e_(k_) « + (k+)-e_(k_))3 



. (23) 



The frequency independent term in eqn. (|23|) corre- 
sponds to the current density i/B introduced in eqn. Q . 
We are interested in linear in q part of density of cur- 
rent. Expanding the integrand up to the first order in 
T^q one can prove by direct calculation that this term 
vanishes. Thus, in the normal state v — as it should 
be in gauge invariant theory (see Section V ancU^) . The 
frequency dependent term determines the current 

e 2 7 2 

jf(u, q) = ieiji— ■^■hq m (hu) 2 Ii m A j (u, q), (24) 



where 



d 3 k 
(2irh) 3 ' 
1 



(£+-£-) 4 \dk m dk 
/a/(e+) ae+ , one-) 



(£+-£-) 3 V dt+ dk n 



<9£_ dk„ 



(25) 



After substitution of the Fourier component of the vec- 
tor potential by the Fourier component of an electric field 
A = cE/iui, we obtain 



jf (w,q) = e lJ ie 2 h 3 ~f 2 u;q m Ii m Ej(u!,q). 



(26) 
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The integral I/ m given by eqn. (|25|) consists of two 
different contributions. The first part of it is determined 
by the difference of the Fermi distribution function for the 
quasiparticles in two bands, another one originates from 
the derivatives of these functions. Performing integration 
over momentum space for the spherical Fermi surfaces in 
the limit 70 &f <C Ef, we obtain 

0/ \ f 1 1 \ e 2 Lo 

j? (w,q) = eyi o - 777 ) _2- 7. g^j(^q) 



24 / ir 2 j k F 



= e 



2 



127r 2 7 ofc F 



^■(w.q). (27) 



Here, in the first line, we mark out the contributions from 
the two parts of the integral (1231) . The corresponding 
gyrotropy conductivity is 



A = i 



2 

e u> 



127r 2 7ofcF 



(28) 



IV. GYROTROPY CONDUCTIVITY IN THE 
SUPERCONDUCTING STATE 

To find the gyrotropy conductivity in the supercon- 
ducting phase with the cubic symmetry, one needs to 
perform summation and integration in the eqn. (1201) us- 
ing G and F in the superconducting phase. The integral 
in eqn. (I2TH) consists of two terms 



d 3 k 
{2irh) : 



■A,[J+_(k,w)-J_ + (k,w)] ) (29) 



where 



m= — 00 

(30) 

and 

00 

J_+(k,u;)=T J2 lG-(K+)G + (K_)-F_(K + )Fi(K_)} 

(31) 

I 



m— — 00 



Transforming the summation into an equivalent con- 
tour integration—, eqn. ([50)) can be written as 



(k, 10) = d> doj' tanh — 

Airi J 2T 

: { [G*iV, k+) - G+(u/, k + )]G A (uj' - u , k_) 
+ [G fl V, k_) - G A (w', k_)]G? (w' + w, k+) 
-[ff( W ',k + )-F4 4 ( w ',k+)]^( w '- w ,k_) 

-[F fl V,k_)-F^( w ',k_)]F^(c' + w ,k + )}. 

(32) 



Here, the Green functions are 



Gf>,k) 
and 

Ff' A ( W ,k) = 
where 



?kJ — ex(k) ± i<5 Tiu; + e A (k) ± iS 

tx(k)A r 1 1 
2e A (k) 



(33) 



+ e\(k) ±18 hu> — e\ (k) ± i<5 

(34) " 



«l(k) 
^(k) 



2^ ± e A (k);' 



(35) 



e 2 A (k)=£ 2 (k) + A 2 , (36) 
and k± = k ± q/2. Taking into account 

Gx-Gi = -27rt[u^(w - ex) + v\8{uj + ex)}, (37) 



A J A 



[d(w-e A )-*(w + eA)], (38) 



after integration with respect to uj', we can rewrite eqn. 
(1321) as: 



tanh 



tanh 



Mk+ 
2T 

e+(k+) 
2T 



tanh 



tanh 



e-(k- 
2T 

e-(k_; 



f A 2 

"24e + (k + )e_(k_) 



, e+(k+) 
tanh +v +; - tanh 
2T 



2T 
e_(k_ 



2T 



/ , e+(k+) , e_ k_ 

tanh -!A_+Z + tanh — 1 '- 

\ 2T 2T 



? 4(k + ) u 2 _(k_) 



« 2 (k + ) v 2 (k_; 



e+(k+) -e-(k-) -w e+(k + )-e_(k_)+ W 
<(k + ) V 2 (k_) <(k + ) U 2 _(k_) 



e+(k+) + e_(k_) -w e+(k+) + e_(k_) + u 
1 1 



e+(k + )-e_(k_)-w e + (k + )-e_(k_)+ W 
1 f 



e+(k+) +e_(k_) - w e+( k +) + e_(k_) + w 



(39) 
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Here we have ignored the shifts in the arguments of the 
phase factors: t\(k± q/2) » t\(k) leading to the small 
corrections of the order of 7ofcp/£F to the main terms. 

For the second term under integral in the eqn. (|29l) we 
have 



J ^(k,w) = (— k, — oj). 

Hence, the integral ([25)) can be rewritten as 
d 3 k 



h = 



rfc;[J+_(k,w) + J + _(k,-w)] 



(40) 



(41) 



It means that we should work with the doubled even 
part of eqn. (|39[) . Expanding the integrand in powers of 
uj after long but straightforward calculations we come to 
the following formula 



As in the normal state, the gyrotropic current consists 
of two different contributions. The first part of it is de- 
termined by the quasiparticles filling up the states in be- 
tween two Fermi surfaces. Another part originates from 
the derivatives of electron distribution functions in two 
bands. The first contribution is not changed in the su- 
perconducting state, at A << 70 But the second 
contribution in the superconducting state is suppressed 
in comparison with its normal value due to the gap in the 
quasiparticle spectrum. Neglecting the terms of the order 
of A 2 /7q/c|, one can substitute under integrals the deriva- 
tives of the particle distribution functions dn(£±) / d£± by 
the derivatives of the quasiparticles distribution function 
df(e±)/de±, where /(e) = (1 -tanh^,)/2. This leads 
to the temperature dependence of gyrotropy coefficient 
A(T) determined by the integral lf m 



d 3 k 
{2nh) 3 
1 



fc ; [n(£+(k+))-n(£-(k_)) 

(M 2 



-2A' 



£+(k+)-£_(k_) 
d 3 k t 



.(k + )-£_(k_))3 

(huj) 2 



A(T) = 
Here Y(T) = f ■ 



87T 2 7 fci: 



IT J cosiJt/2T ^ S ^ aB Yosida function. It 
is useful to compare this formula with the corresponding 
normal state equation (|27|) . 

Besides of this term there is another part of the current 
originating of the integral 1} , It yields the coefficient 



(1 - Iy(t)). 



(48) 



(2nh) 3 l (£+(k+) 
/tanhS±^ 

X ; 



2T 



+ £_(k_))(£+(k + ) 

tanh^^ 



C-(k-)) 3 



2T 



where 



n(0 



1 



e_(k_) 



— tanh — 
e 2T 



(42) 



v{T) 



e 2 7ofcf 
6n 2 h 2 i 



(1 - Y(T)) 



(49) 



(43) 



such that the total gyrotropy coefficient A(T) in the su- 
perconducting state is 



is the distribution function of electrons over energies. 
At A/'yokp <C 1 the second integral is obviously much 
smaller than the first one. So, we come to the expression 



A(T) 



270 h 



1 - \y{t) 



d 3 k 
{2nh) 3 
1 



«+(k+)) 



n(f_(k_))] 
(hu) 2 



e + (k + )-e-(k-) (e + (k + )-e-(k-)) 3 



(44) 



which has the same form as the corresponding formula 
for the normal state (|2"3")l . 

Expanding the integrand up to the first order in — ' • 
we obtain for the current given by eqn. (HI 



dk 



e 2 7o 2 



jf (w, q) = ieijl —^hq m [J/ ro + (Huj) lf m ]Aj ( U , q), (45) 



(l-Y(T)) 
.(50) 

To find a relationship of the gyrotropy conductivity 
with observable optical properties one has to develop 
electrodynamic theory of noncentrosymmetric metals. 



V. OPTICAL PROPERTIES OF 
NONCENTROSYMMETRIC METAL 

A. Dispersion law 

To derive the light dispersion law we start from the 
Maxwell equations 



I lm — 



d 3 k 
{2nh) 3 
1 



»(£+)- n(£_) (di 



(£f- 



£-) 2 



dk„ 



+ 



£+ - £- V <9£+ dk n 



<9£_ dk n 



1 lm — 



+ 



d 3 k 
{2TTh) 3 ' 
1 

(C+-C-) 3 V dt+ dk n 



,n(f + )-n(e_) 



(9fc„ 



rotB 



rotE 



4tt. 

— J , 

c 

1 dB 



(46) 



supplied by the density of current expression 



+ 



\dk„ 



dk„ 



e <9E 



crE + ArotE. 



(51) 
(52) 

(53) 



<9£_ dkn 



.(47) 



The first term here corresponds to the dispersionless part 
of the displacement current. The second one is the con- 
ductivity current written at infrared frequency region 
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u) > vp/S, lot > 1, where the current is locally related 
with an electric field, 5 is the skin penetration depth. 
The last one is the gyrotropy current 



j g = ArotE. 



(54) 



As before we discuss the metal with the cubic point sym- 
metry. 

Eliminating the magnetic induction, we obtain 



V 2 E 



e <9 2 E AnadE 4ttA SrotE 



c 2 dt 2 c 2 dt c 2 dt 
Taking solution for the circularly polarized wave 

E = (x ± iy)E Q e i( - kr - ut ^ 

we come to the dispersion relation 

• 2 Airiaui AitiXujk 



k A = 



ecu 



± 



(55) 



(56) 



(57) 



It is worth to be noted that for a media with time reversal 
breaking one has to substitute here a — > <j± = <j xx ±ia xy , 
where a xy is the Hall conductivity. 

In neglect the gyrotropy term the complex index of 
refraction 

ck 

N = — = n + in 

is expressed through the diagonal part of complex con- 
ductivity a = a' + ia" by means of the usual relations 



Ana 1, 



k = e — 



2nn 



4na' 



The gyrotropy term leads to the difference in the refrac- 
tion indices of clock wise and counter clock wise polarized 
light. In the first order in respect to A = A' + i\" the 
refraction index is 



AT* = n + in ± 



2iri\ 



(58) 



Hence, the differences in the real and imaginary parts of 
the refraction indices of circularly polarized lights with 
the opposite polarization are 



An = n + — n_ = — 



4ttA'' 



An = 



4ttA' 



(59) 



(60) 



In the superconducting state the gyrotropy current 
([54")) has more general form given by eqn. Hence, 
we should use the more general formula for the current 



e <9E 



crE + ArotE + isB. 



(61) 



Then repeating all the calculations we come to the same 
results ([57l) - ([6"0]) modified by the substitution 



A 



A 



A lCV 

LO 



(62) 



We remind that the superconducting state current den- 
sity given by eqn. (|61[) is worth to use at the high 
frequencies co > vf/S where the inequality Tiuj >> A 
is certainly valid. Here, A is superconducting energy 
gap. In the low frequency limit TkJ < A one should also 
take into account the London density of current jl = 
— (c/47r5|)(A — hcVtp/2e). The interplay between the 
London current jz, and the Drude current = ct(u;)E is 
discussed in the textbook^ 



B. Magnetic moment 



The magnetization in gyrotropic media is 



m = 1ae, 



(63) 



such that rotation of the magnetization is equal to one 
half of gyrotropy part of the current density 



1 



c rotM. 



(64) 



The relationship between the density of gyrotropy cur- 
rent and the magnetization is a general property of non- 
centrosymmetric materials (see also^). Both of them can 
be obtained from the gyrotropy term in action 



1 



^d 3 r(AE)B. 



By variation of action in respect of — B and — A/c, taking 
into account that A is an odd function of derivative d /dt 
and making use the definitions E = —{l/c)dA/dt and 
B = rot A, we come to M given by eqn. (1531 and j g given 
by eqn. (|54p correspondingly. 

All these considerations are valid both for the normal 
and as well for the superconducting state. However, in 
the latter case the gyrotropy action 



1 

2c 



dtd 3 r |(AE)B 



A 



he. 

Ye 



B 



(65) 

contains one extra term which is absent in the normal 
state due to the gauge invariance. The corresponding 
expressions for the magnetic moment and gyrotropy cur- 
rent are 



M 



1 

2c' 



-AE 



1 

— i 

2c 



Tic. 
2e" 



A rotE + vB. 



(66) 



(67) 
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C. Boundary conditions 



D. Reflection coefficient and the Kerr effect 



To consider the problem of light reflection normally 
incident to the flat surface of noncentrosymmetric metal 
we need to find the relations between the wave amplitude 
propagating inside (z > 0) the material 



£■» _ jj _iw(JVa/c-t) 



(68) 



and the amplitudes of incident and reflected waves out- 
side it 



■gout _ g iu(z/c-t) _|_ g -to<(*/ 



c+t) 



(69) 



We have 



that is 



Ein Tj\out 



Eo = Ei + E2. 



(70) 



At the same time from the difference of the Maxwell equa- 
tions (|52[) inside and outside of material we obtain 

(rotE m - rotE out ) z=0 = --^(B m - B OMt ) 2=0 . (71) 

c at 

The difference of the magnetic inductions at the bound- 
ary is given by the jump of magnetzation 



(B ra - B OMt ) 2=0 = 4^M 2=0 . 



(72) 



In the stationary magnetic field parallel to the surface 
of the metal H ou * = H x x this equation yields 

B l x nt {z = 0) = H x , B y nt (z = 0) = —A y nt (z = 0). 

In the normal state where k = the boundary condi- 
tions add nothing special to the centrosymmetric case. 
In the superconducting state the solution of the London 
equations supplied by these boundary conditions results 
in quite unusual helical field distribution found in the 
paper »ii 

For the light incident to the metallic surface using (I7T1) , 
([72]), and dSU) we obtain 



4tt /. d 



(rotE m - rotE otlt ) 2=0 = ( ^ " <* ) K=o (74) 

Substituting here eqns. (jB"5|) . we come to 

z x (iVE - Ei + E 2 ) = — AE (75) 

c 

For the combinations E^ = E x ± iE y of the electric field 
component this relation can be rewritten as 



iV ± ± Eq -Ef + Ef = Q (76) 



The equations (|70|) and ([76]) allow express the ampli- 
tudes of reflected wave through the amplitude of the in- 
cident wave. We have for reflection coefficient 



1 - N± =F 



2niA 



Ef 1 + 7V± ± 
where the refraction index is 

2?riA 



2ttiA ' 



(77) 



N- 



± 



■ in ± 



(78) 



Now one can rewrite eqn. (1771) in more habitual form 

1 - A>± 



R ± = 



1 + JV± 

where an effective refraction index is 

47TlA 

Af± = n + zk ± , 



(79) 



(80) 



and the effective differences in the real and imaginary 
parts of the refraction indices of circularly polarized lights 
with the opposite polarization are 



An = 



Ak = k + — k_ = 



8ttA" 



8ttA' 



(81) 



(82) 



Making use these definitions we can apply the standard 
procedure^ to calculate the Kerr rotation for linearly po- 
larized light normally incident from vacuum to the flat 
boundary of a medium. The light is reflected as ellip- 
tically polarized with the major axis rotated relative to 
the incident polarization by an amount 



e 



(i 



n 2 )Ah + 2n,KAh 



(l-n 2 + K 2 ) 2 + (2n K ) 2 



(83) 



VI. THE KERR ROTATION 

To find the Kerr rotation in the normal state let us sub- 
stitute the eqn. (gSJ) in eqns. (JHTJ), ([52)1. We find Ak = 
and An expresses through ratio of the light frequency to 
the band splitting 270 &f as 



An = - 



2a Tiuj 
3ir 7 fc F 



(84) 



Here, a = e 2 /he is the fine structure constant. 

We limit ourselves by the frequencies not exceeding the 
band splitting jokp. Although the latter is not known 
for many noncentrosymmetric materials, one can expect 
it is about hundred Kelvin or in the frequency units 



8 



~ 10 13 r ad I sec As an example we consider the sit- 
uation when the frequency of light is of the order of this 
value and larger than the quasiparticles scat tering rate 
(clean limit): 1 << ujt < w p r, where u p = \J Aime 2 Jm* 
is the plasma frequency. In this frequency region the real 
and imaginary part of conductivity are a' « uj 2 /4ituj 2 t 
and a" m uj 2 /4ituj. Then, one can find 2tik ~ oj 2 /oj 3 t 
and n 2 — n 2 « ui 2 /ui 2 . Thus, for the Kerr angle we ob- 
tain 



2a hu 2 

37T J k F UJ 2 T' 



(85) 



So, the Kerr angle in noncentrosymmetric metals can 
have measurable magnitude, in particular if we compare 
it with the Kerr angle of the order of 6 x 10~ 8 rad mea- 
sured in the superconducting S^RuO^ by the Stanford 
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group.— 

For An in the superconducting state we obtain 



2a 

An= 

7T 



TlUJ 

2j k t 



(86) 

Finally, for the Kerr angle in the same frequency interval 
as for the normal state we have 



and wave vector in noncentrosymmetric metal. The ob- 
tained general formula valid both in the normal and in 
the superconducting state was applied to the calculation 
of observable physical properties in the frequency interval 
smaller than the band splitting Huj < jokp. The latter 
in its turn was supposed to be smaller than the Fermi 
energy 70 &f < £f- The calculations was performed in 
the clean case lot > 1, that is, in particular, important 
to neglect the vortex corrections. We did not discuss 
the anomalous skin effect assuming that the wave length 
does not exceed the skin penetration depth S > vf/lj. 
In the normal state the current contains the gyrotropic 
part which is odd function of the wave vector and the fre- 
quency. It presents a sort of displacement current orig- 
inating of band splitting in noncentrosymmetric metal. 
In the superconducting state there is an additional part 
of the gyrotropy current proportional to magnetic field. 
The temperature dependence of gyrotropy conductivity 
in the superconducting state was found. As an example 
the Kerr rotation for the polarized light reflected from the 
surface of noncentrosymmetric metal with cubic symme- 
try is calculated. 
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